Abstract We examined the sky confusion noise in 40 sky 
regions by analysing 175 far-infrared (90-200 /Ltm) maps 
obtained with ISOPHOT, the photometer on-board the 
Infrared Space Observatory. For cirrus fields with (B) > 
5MJysr _1 the formula based on IRAS data (Helou & Be- 
ichman, ( 1990] )) predicts confusion noise values within a 
factor of 2 to our measurements. The dependence of the 
sky confusion noise on the surface brightness was deter- 
mined for the wavelength range 90 < A < 200 /im. We 
verified that the confusion noise scales as Noc (B) 15 , in- 
i— ( ] dependent of the wavelength and confirmed N oc A 2 5 for 

■ A > 100 /an. The scaling of the noise value at different 
| separations between target and reference positions was 

■ investigated for the first time, providing a practical for- 
| mula. Since our results confirm the applicability of the 

■ Helou & Beichman ( 1990 ) formula, the cirrus confusion 
j noise predictions made for future space missions with tele- 

■ scopes of a similar size can be trusted. At 90 and 170 /j,m 
I | a noise term with a Poissonian spatial distribution was de- 

■ tected in the faintest fields ((B) < 3 — SMJysr -1 ), which 
C<") ' we interpret as fluctuations in the Cosmic Far-Infrared 
■"xt" Background (CFIRB). Applying ratios of the fluctuation 
' , _ , ' | amplitude to the absolute level of 10% and 7% at 90 

I ■ and 170 ^m, respectively, as supported by model calcu- 
I ' lations, we achieved a new simultaneous determination of 
f^*) • the fluctuation amplitudes and the surface brightness of 
"^h \ the CFIRB. The fluctuation amplitudes are 7±2mJy and 
O ■ 15±4mJy at 90 and 170 /j,m, respectively. We obtained 
J^;a CFIRB surface brightness of B = 0.8±0.2 MJysr" 1 
(id„ = 14±3 nWm sr ) at 170 /jm and an upper limit 
^ of 1.1 MJysr" 1 (v\ v = 37nWn:r 2 sr- 1 ) at 90/zm. 

: 
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1. Introduction 

Sky confusion noise causes an uncertainty in the determi- 
nation of the source flux, due to the variation of the sky 
brightness between the target (on-source) and reference 
(off-source) positions. This noise cannot be overcome by 
longer integration times, thus it constitutes a basic limi- 
tation for the detection of very faint point sources. 

At far-infrared wavelengths the two major components 
responsible for the sky confusion noise are dust emission 
from irregularly sh aped interstellar clouds, the "galactic 
cirrus" (Low et al., 1984 ), and a component of the Extra- 
galactic Background built up from the accumulated light 
of faint unresolved galaxies along the line of sight (Guider- 
doni et al., |1997[ ). 

The confusion noise due to cirrus can be characterized 
by a formalism applied first by Gautier et al. ( |1992|) for 
the IRAS 100 /im scans. They computed the second order 
structure function for a far-infrared map : 
B(x-0)+B(x + l9) 



S(0) = B(x) - 



(1) 



where B is the sky brightness, x is the location of the tar- 
get, 9 is the separation between the target and reference 
sky positions and the average is taken in spatial coordi- 
nates over the whole map. The noise due to sky brightness 
fluctuations, N, is defined as: 

N(0) = y/^tffj x (2) 
where fl is the solid angle of the measuring aperture. 

The measured confusion noise depends both on the size 
of the measuring aperture D and on the angular separation 
9 between the on-source and off-source positions. However, 
if the sky brightness distribution does not show any char- 
acteristic scale length (e.g. random fluctuations or fractal 
structure), the only important factor is the 9/D ratio. For 



their study of the IRAS 100 /im maps Gautier et al. ( 1992 ) 
defined a standard measurement configuration where the 
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on-source position is bracketed by two off-source positions 
separated by ± 9 (Fig. Q) and the ratio is fixed to the res- 
olution (Nyquist) limit 6> = # m i n = 2D. 
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Figure 1. Measurement configuration to compute the 
structure noise. Gautier et al. ( 1992j ) and Helou & Beich- 
man (1990) derived their results for the resolution limit 
2D. In Sect. 3.2 we also consider the cases of 



/ min — *xj. In Sect. 
9 = k-D,k = 2...7 



Gautier et al. (1992) pointed out that the Fourier 



power spectra of the IRAS 100/im scans were generally 
well represented by a power law, and that the structure 
noise can be linked to the parameters of the power spectra 
by the following relation: 

Noc^xpj (3) 

where d^ 1 is the spatial frequency, d^ 1 is a reference spa- 
tial frequency, Po is the Fourier power at d^" 1 , and a is the 
spectral index of the Fourier power spectrum. It was also 
found that most scans can be described by the empirical 
relationships a«3 and Po cx (B) 3 . 



Based on the IRAS results Helou & Beichman (1990) 



(H&B) proposed a practical formula to predict the cirrus 
confusion noise, assuming that the empirical relationships 



of Gautier et al. (1992) are valid at other wavelengths, too. 



They also took into account that at the resolution limit 
of the telescope d/do can be expressed by the resolution 
parameter of Fraunhofer diffraction, A/D t , where A is the 
wavelength of the measurement and D t is the diameter 



Cs. Kiss et al.: Sky Confusion Noise in the Far-Infrared 



3 



of the telescope primary mirror. Under these terms they 

expressed the cirrus confusion noise as: 

n h & b / A y.*,p t y2.*, (Ba ) y-5 

lmJy ' V 100 Aim/ Vim/ VlMJysr- 1 / 1 ' 
Here B\ is the surface brightness at the wavelength of 
the observation. The formula shows that the confusion 
noise depends on the wavelength in two separate ways: 
(1) via the variation of the surface brightness (spectral 
energy distribution of the emitting medium) and (2) via 
the dependence of the resolution parameter on A. This re- 
lation which is valid for the standard measurement con- 
figuration as defined in Fig. |l| has been widely used for 
preparing far-infrared space telescope observations, and 
applied even for other configurations. The first noise anal- 
ysis of four I SOPH OT maps was performed by Herbst- 
meier et al. ( 1998 ), carrying out the first high spatial 
resolution study of the galactic cirrus at A>100/im. They 
proved that the ISOPHOT sensitivity at far-infrared wave- 
lengths is limited by sky confusion noise rather than by 
instrumental uncertainties. 

The confusion noise caused by the CFIRB has dif- 
ferent properties from the cirrus confusion noise. It can 
be represented by a Poissonian spatial distribution, there- 
fore the confusion noise value is independent of the sep- 
aration between the target and the reference positions. 
Ackermann et al. (1992) investigated the sensitivity lim- 
its, including the effect of the CFIRB, for ISOPHOT 
and predicted that the long wavelength observations 
would be limited by galaxy confusion in the regions of 
faintest cirrus. Recently Lagache & Puget (200C), Mat- 
suhara et al. ( 2000|) an d Juvela et al. (for an overview 
see Lemke et al. T POOOD ) analysed ISOPHOT maps and 
separated the cirrus and extragalactic components. They 
announced the detection of the CFIRB with a fluc- 
tuation power in the range of 5 — 12xl0 3 Jy 2 sr _1 at 
170 and 1 80 /im , close to the predictions of Guider- 
doni et al. ( 1997] ). This leads to a galaxy confusion noise 
limit of 13-22 mJy for ISOPHOT at this wavelength. 

For this paper we analysed 175 suitable far-infrared 
maps from the ISO Archive. Our goals w ere ( 1) to test 
the applicability of the Helou & Beichman ( 1990| ) formula, 
especially at A > 100 //m; (2) to determine the relative 
contributions of the galactic cirrus and the extragalactic 
components in the faintest regions of the sky; (3) to de- 
rive an easy-to-use formalism for predicting the total sky 
confusion noise for instruments on ISO, SIRTF and HER- 
SCHEL and for various measurement configurations, and 
(4) to determine the Cosmic Far-Infared Background Ra- 
diation via its fluctuations. 

2. Observations and data analysis 

2.1. Selection of ISOPHOT maps 

We s elected 175 maps from the ISO Archive (Kessler et al., 
2oobT ), observed in the 90 to 200 /im wavelength range, cov- 
ering 40 different sky regions. All maps were obtained with 



the ISOPHOT instrument (Lemke et al., 1996) on-board 
the ISO satellite (Kessler et al., |1996|) in the PHT22 raster 
observing mode (Laureijs et al., 2000| ). Our selection cri- 
teria were that the maps should be larger than 5x5 raster 
points, corresponding to ~ 8' x 8' for the C100 (3x3 pixel 
array, 46" pixel size) and C200 (2x2 pixel array, 92" pixel 
size) detectors. Some maps were observed with full over- 
sampling, i.e. the field was redundantly covered without 
gaps between detector pixels. The final map images are a 
composite of the individual array images on the different 
raster positions. The size of the image pixels is always a 
detector pixel. 

We excluded maps with obvious individual structures 
(resolved stars, galaxies, planetary nebulae, etc.). The 
number of selected maps is 94, 4, 65 and 12 using the 
C1.90, C1JL00, C2.170 and C2.200 filters, respectively, 
and they cover a wide range of surface brightnesses (~ 1- 
100 MJysr -1 ). A more detailed description of the maps 
will be presented in a forthcoming paper (Kiss et al. 2001, 
in prep.) dealing with the spatial structure of the extended 
emission. 

2.2. Data reduction 

Basic data reduction was performed using the ISOPHOT 
Inter active Analysis software, PIA V8.2 (Gabriel et al., 
1997 ) with the standard batch mode set-up. Flat-fielding 
had to be applied, which corrects for the residual respon- 
sivity differences of the individual detector pixels. We 
chose the First Quartile Normalization method, which 
uses the first quartiles of the maps' brightness distribution, 
computed for each detector pixel, for normalization. The 
reliability of this statistical flat-field method was tested 
in the case of oversampled maps, were redundant obser- 
vations of the same sky position by each detector pixel 
are available. These tests confirmed that the First Quar- 
tile Normalization gives accurate flat-field values, and was 
therefore adopted for this analysis. 

The maps still contain the zodiacal light (ZL) emis- 
sion, and a contribution from the CFRIB. The zodiacal 
emission does not contribute to the sky confusion noise, 
because its distribution is sm ooth on the scale of our 
maps (Abraham et al., 1997). Therefore, we can elimi- 
nate it from the total brightness. First, we determined 
the non color corrected ratio of the ZL-to-total emis- 
sion at the DIRBE wavelengths, using the COBE/DIRBE 
weekly maps and the COBE/DIRBE sky and zodi atlas 
(Hauser et al. Ql998| ), Kelsall et al. fll998| )). We colour cor- 
rected the ZL component adopting a spectral energy dis- 
tribution (SED) shape of a black body with a temperature 
of 270 K (Leinert et al. ( pOOl) ) find 255 K < T ZL < 297 K) . 
The difference between the uncorrected total and ZL emis- 
sion was color corrected by assigning a cirrus SED to 
it, assuming a 20 K modified black body with v 2 emis- 
sivity law. Then interpolating the corrected ratios in be- 
tween the COBE/DIRBE wavelengths, we estimated the 
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ratios of the ZL-to-total emission for the ISOPHOT filter 
central wavelengths. The absolute ZL contribution in the 
ISOPHOT maps was computed by multiplying this ratio 
with the average total brightness measured by ISOPHOT 
(see also Heraudeau et al., 200l[ ). 

The remaining surface brightnesses include the CFIRB 
and the cirrus emission. In Sect. || we calculate the sky 
confusion noise and in Sect. |^2] we separate the cirrus 
and Extragalactic Background fluctuations, and derive the 
value of the CFIRB. 



2.3. Instrument noise 

The determination of the instrument noise is crucial in 
order to get the real value of the sky confusion noise. In 
addition to the basic noise components (read-out noise, 
dark current variations, etc.), we also include in the in- 
strument noise uncertainties related to the generation of 
the final maps (flat-field). 

There are several ways to determine the instrument 
noise. The fundamental difference between them is the 
time scale over which the noise estimation was taken. We 
examined the following four methods: 

— Ramp-noise: Each measurement is composed of many 
individual integrations on each raster position. A set 
of non-destructive read-outs builds up an integration 
ramp whose slope provides the signal of this integra- 
tion (Lemke et al., 199£). From each measurement we 
created two maps using the even and odd integration 
ramps separately. The instrument noise was calculated 
from the difference between the 'even' and 'odd' maps. 
This noise estimate provides information on the detec- 
tor stability on a time scale of seconds, the typical 
difference between adjacent ramps. 

— PIA-noisc: PIA provides a determination of the signal 
uncertainty of the individual raster positions, which 
can also be adopted as an estimate of Ni nst . It reflects 
the error propagation of random error components in 
(1) the linear fitting of the integration ramps, and (2) 
the averaging of all signals per raster point [j. The time 
scale here is the time spent on an individual raster 
position, typically in the order of a minute. 

— Flat-field-noisc: This noise value was calculated from 
the variation of the brightness at the same sky posi- 
tion in maps created from individual detector pixels 
(i.e. in case of pixel redundancy), after correction by 
the appropriate flat-fielding factors. Since these fac- 
tors are kept constant for the whole map, although 
the detector pixels show different time dependent tran- 
sient curves, this instrument noise estimate includes 
the effects of the long term changes in the detector be- 
haviour. Therefore the time scale is the observational 
time of the whole map, typically in the order of an 



hour. Since the transient behaviour depends on the 
illumination level, the flat-ficld-noise may also depend 
on the brightness of the map. 
— Repetition- noise: If a sky region has been observed sev- 
eral times, then it is possible to derive a noise value 
from the comparison of the individual maps. The typ- 
ical time scale here is the difference in the observation 
dates. It provides a measure of the measurement re- 
producibility. 

The ramp- and PIA-noise could be calculated for 
each map, but the flat-field-noise could be evaluated only 
in maps with full oversampling (see Sect. 2.1) and the 
repetition-noise was only computed for one particular sky 
region (Marano 1 field) consisting of several submaps, each 
observed four times. A comparison of the different noise 
estimates showed that the ramp-noise provided the lowest 
Nmst values, typically a factor of 2 lower than the PIA- 
noise. The flat-field noise was higher than the PIA one, 
and the repetition noise was found to be nearly identical 
to the flat-field-noise. Since in the case of the structure 
noise calculation (Eq. [I]) the target and reference posi- 
tions were observed at times typically separated by several 
minutes, the correct instrument noise should lie between 
the PIA- and the flat-field-noise values. To be on the safe 
side, we chose the value of the flat-field-noise, although it 
may represent a somewhat conservative estimate, which 
might cause a slight underestimation in the final confu- 
sion noise values. Since this noise value is not available 
for each map, we took the uncertainties provided by PIA 
and scaled them according to the mean ratio of 1.35 and 
1.65 of the flat-field- and PIA-noise values for the C100 
and C200 detectors, respectively. 

The typical instrument noise values we found (see 
Fig. ||) are somewhat higher than that of other authors, 
e.g. a factor of 2 compared to Lagache et al. (200C). This 
discrepancy can be traced back to the co-addition of four 
independent maps in their case, reducing the final noise 
values by this factor. 



The on-line PIA-manual is available at the URL: 



" |http:/ /www. iso.vilspa.esa.es/manuals/PHT/pia/pia.html 



Following Herbstmeier et al. ( 1998 ) we assumed that 
the sky brightness fluctuations and instrument noise con- 
tributions are statistically independent, therefore the mea- 
sured structure noise can be expressed as: 

N mcas < N 2 + 2 • N 2 „ st (5) 

where N moas is the measured structure noise, N is the real 
structure noise and Ni ns t is the instrument noise. We used 
this formula for the subtraction of the instrument noise 
from the measured noise, by assuming equality between 
the two sides of the expression. In order to test this as- 
sumption we compared the confusion noise values of four 
C2_170 maps, where the same sky region was mapped re- 
peatedly. The confusion noise must be the same and only 
the instrument noise values can differ in repeated obser- 
vations. Assuming equality in Eq. |5| we obtained identical 
confusion noise values for all four measurements and for 
any separations. 
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Figure 2. Instrument noise values measured on the 
ISOPHOT maps. In the case of oversampled maps the 
flat-field noise, otherwise the scaled up PIA-noise was 
adopted. The mean scaling factor between the PIA- and 
the flat-field noise values are 1.35 and 1.65 for the C_100 
and C_200 detectors, respectively, (a) C100 detector: 
Cl_90 and Cl_100 filters (plus signs and diamonds, re- 
spectively), (b) C200 detector: C2.170 and C2.200 filters 
(black dots and crosses, respectively) 



3. Results 

3.1. Confusion noise at the resolution limit 

Fig. H presents our results on the sky confusion noise (af- 
ter removing the instrument noise using Eq. |^) measured 
with the four ISOPHOT filters COO, C1.100, C2.170 and 
C2_200. We us ed the convention 9 = 2D, defined by Gau- 
ticr et al. ( |1992| ) (see Sect. Q). Following the H&B formula 
(Eq. ||) we fitted the data points for each filter separately, 
assuming a power-law relationship between the confusion 
noise and the average brightness of the field, but also al- 
lowing a constant factor: 

N(0 mi „) / (B) ^ 

— i. ii + Lax^- 



lmJy ~ u ' ~ x " VI MJysr- 1 ) ^ 
The coefficients Co, Ci and r\ are listed in Table [I]. In 
the case of the Cl_100 and C2_200 filters no appropriate 
determination of the Co offset was possible due to the lack 
of faint regions. Eq. ^ can be used to predict the confusion 
noise for the ISOPHOT filters studied here. Our confusion 



Figure 3. Confusion noise at the resolution limit # m i n vs. 
surface brightness of the fields, (a) ISOPHOT C100 filters, 
C1.90 and C1.100. (b) ISOPHOT C200 filters, C2_170 and 
C2_200. The symbols are the same as in Fig. ^ 



filter 


C1.90 


C1.100 


C2.170 


C2.200 


Co 


6.6±1.9 




10.2±2.5 




Ci 


0.7±0.3 


0.6±0.1 


3.0±0.3 


4.4±0.2 


V 


1.55±0.12 


1.53±0.28 


1.47±0.11 


1.56±0.19 


Pmin 


92" 


92" 


184" 


184" 


Q (10~ 7 sr) 


0.6469 


0.7030 


2.6438 


2.8120 



Table 1. Results of the fits using the measured values 
of the sky confusion noise at the resolution limit N(6* m ; n ) 
and the average surface brightness of the field (B) (see 
also Eq. ||and the text for details). Separation limit and 
effective solid angle values are also presented. In the case 
of the Cl_100 and C2_200 filters no appropriate determi- 
nation of the Co offset was possible due to the lack of faint 
regions. 



noise values measured at 170 /mi are in agreement with 
a value of 45mJy obtained by Dole et al. (2001) in the 



FIRBACK regions, taking into account that their value is 
a source flux confusion limit, which differs from ours by 
the footprint-fraction and is therefore a factor of 2 higher. 
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Figure 4. Confusion noise at the resolution limit deter- 
mined by this stud y as compared to the prediction of 
Helou & Beichman ( |l99d| ). The symbols are the same as 
in Fig B 



3. 2. Comparison with the H&B formula 

In Fig. U we compare our results with the predictions of 
the H&B formula (Eq. [|). In the surface brightness range 
where the cirrus emission dominates (5-30MJysr _1 ) the 
H&B formula predicts confusion noise values within a fac- 
tor of 2 of our results. The scatter cannot be attributed to 
measurement uncertainties, and probably it reflects the 
physical differences among the fields. For the brightest 
fields the H&B formula seems to systematically overesti- 
mate the measured confusion noise values. Such regions, 
however, contain molecular clouds, therefore they may 
have different spatial structure from a cirrus region and 
should be described by a different power law. The figure 
also shows a large discrepancy at low surface brightness 
which could well be treated by the introduction of an off- 
set (see above). A possible explanation of the deviations 
in the faint range is discussed in Sect. 4.2. 



Although the H&B formula gives rather good predic- 
tions for the surface brightness range it was developed for 
(5-30MJysr _1 ), its coefficients, especially the three expo- 
nents, rely partly on the analysis of IRAS 100/im data, 
partly on theoretical considerations. Our data set offers 
a unique chance to test some of these exponents for the 
first time. The dependence of the confusion noise on the 
surface brightness was described by H&B using a power 
law with an exponent of 1.5 (Eq. ||), independently of the 
wavelength. This exponent rj was determined by our fits 
for four wavelengths in Sect. |3.1[ The measured values in 
Table |l| are very close to the value of 1.5, proving the 
wavelength independence of this exponent. 



5.0 



2.0 



1 .0 



0.5 



100 200 
wavelength [jitm] 



Figure 5. Ci coefficients of the fits in Table |l| vs. the 
appropriate wavelength (logarithmic- logarithmic plot). 



In the next step we tested the dependence of the confu- 
sion noise on the resolution parameter A/D t , also assumed 
to be in the form of a power law with an exponent of 2.5 
by H&B (Eq. |). Since the diameter of the telescope mir- 
ror D t is the same regardless of the wavelength, in Fig. || 
we plotted the wavelength dependence of the Ci coeffi- 
cient from Table g. The fi gure confirms the assumption 
of a power law, and gives an exponent of 2.53±0.31, in 
agreement with the value in Eq. g. The dependence of the 
confusion noise on the diameter of the telescope (Eq. 
was not possible to check, since the telescope mirror of 
ISO had the same size as that of IRAS. Our results ver- 
ify exponents of the H&B formula for the first time by 
observations over a wide range. 

3.3. Confusion noise for larger separations 

The H&B formula provides information only on the mea- 
surement configuration presented in Fig. |l|, at the resolu- 
tion limit 9 m i n . However, in many applications the beam 
separation may be different. In most ISOPHOT measure- 
ments the observer was allowed to specify any separation 
up to 3'. Staring on-off measurements and sp arse m aps al- 
lowed even larger separations (Laureijs et al., |2000|) . Prop- 
erties of the sky confusion noise for these cases were not 
investigated so far. 

As indicated in Sect. [I], the ratio of the confusion noise 
values measured at different separations is sensitive to the 
spatial distribution of the dominant component of the 
emission (galactic cirrus or distant galaxies). We found 
that the dependence of the sky confusion noise on the 
separation between target and reference sky positions can 
be described by a simple power law: 

N(q-# min ) = N(0 min )xq^ (7) 
where q=l,1^...3^ and 7 is constant for a specific map. We 
fitted 7 for all maps. The results are presented in Fig. |[ 
which suggests that the 7 value is a filter-independent pa- 
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Figure 6. 7 values of the N(q • ^min) = N(6> m i n )-q 7 fits for 
four ISOPHOT filters versus the surface brightness of the 
fields (see the text for details). The symbols are the same 
as in Fig. || 

rameter of the field, depending only on the surface bright- 
ness. For faint fields 7 is close to zero, while it increases 
strongly with increasing surface brightness. We found that 
the behaviour can be approximated by the following rela- 
tion (dotted curve in Fig. ||): 

_ J 0.02 if (B) < 3 MJysr" 1 , . 

7 ~ i 0.65xlogi (B) - 0.26 if (B) > 3 MJysr" 1 (NI 
Eq. and || together with Eq. || and the coefficients in 
Table [l] provide a practical tool for the estimation of the 
confusion noise for a large variety of measurements. 

4. Discussion 

Variations 0/7 

As shown in Sect. |3.3| , for low surface brightness fields 
the 7 values are close to zero. In those regions the noise 
distribution differs from that expected for the galactic cir- 
rus (see e.g. Ackermann et al., 1992)), originating from its 
multifractal structure. 7 « is typical for a Poissonian 
distribution. 

Higher 7 values in brighter regions may be caused by 
the cirrus structure. Despite this general trend the scat- 
ter of 7 at higher surface brightness is relatively large and 
there are a few regions measured by the CI_90 filter with 
high 7 values at low brightness. These discrepancies are 
probably caused by real physical differences between the 
fields. The 7 parameter depends on the spectral index of 
the Fourier power spectrum of the field (cf . Eq. ||) . This 
relationship will be analysed in detail based on experimen- 
tal data in a forthcoming paper (Kiss et al. 2001, in prep.). 
Differences in physical properties may originate in differ- 



ent chemical/dust composition, spatial structure, gas-to- 
dust or molecular-to-neutral gas ratio. 

4-2. Fluctuations due to the Extragalactic Background 
4.2.1. Findings from this work 

When fitting the sky confusion noise as a function of the 
surface brightness in Sect. 3.1 a constant term Co was 



allowed for. At both 90 and 170 /im definite positive val- 
ues were obtained at the 3-4<r significance level. Assum- 
ing that the cirrus confusion noise follows the power law 
behaviour also at very low surface brightness, this offset 
cannot be attributed to cirrus. The spatial distribution of 
the confusion noise below ~3 MJysr -1 has different char- 



acteristics than in brighter cirrus fields (see Sect. 3.3 and 
Fig. |6|). The likely origin of this noise component, as sug- 
gested by its invariability with the surface brightness and 
its Poissonian spatial distribution, would be the fluctua- 
tion due to the Cosmic Far- Infrared Background (CFIRB). 

If this interpretation is correct, than Eq. |^ provides a 
new method to determine the value of the CFIRB fluctu- 
ation. The features of this method are (f) the subtraction 
of the cirrus component via its dependence on the sur- 
face brightness and (2) that this dependence is calibrated 
simultaneously on the brighter cirrus fields of the same 
database. An advantage of this method is the utilization 
of the largest ISOPHOT database for the determination 
of the CFIRB fluctuations. The results are not sensitive 
to characteristics in individual fields. On the other hand, 
any uncertainty in the surface brightness calibration of the 
cirrus would introduce an uncertainty in the fluctuation 
value, too. 



The values given by the fitting procedure in Sect. 3 
for the Co coefficients are 6.6±1.9 and I0.2±2.5mJy at 90 
and 170 fim, respectively. In this fitting the surface bright- 
ness values still included a contribution from the CFIRB 
itself. Subtraction of the CFIRB increases the fluctua- 
tion amplitudes. Varying the CFIRB surface brightness 
Bo within the assumed range of O.i to i. 5 MJysr -1 and 
subtracting it from the total surface brightness we used 
Eq. ^ to derive the corresponding Co coefficients which 
were then transferred to surface brightness fluctuations 
<5B as <5B = Co/$l\, where f2> is the effective solid angle 
of a detector pixel (see Table [l]). 

The results are shown in Fig. [?]. As it can be seen in 
this figure, the dependence of 5B on the assumed value of 
Bo is not very strong, therefore, regardless of the real value 
of the CFIRB, a well confined value of SBo can be derived. 
We obtained C = 7±2mJy (SB Q = 0.11±0.03 MJysr" 1 ) at 
90/mi and C = 15±4mJy (SB =0.06 ±0.02 MJysr" 1 ) at 
170/im. 

In most models CFIRB fluctuations are caused by 
galaxy clustering, and, due to the flat power spectrum, 
the 5Bo/Bo value is con stant in the 1° to 5' res olutio n 
range (Kashlinsky et al. ( |l996[) , Haiman & Knox ( |2000D ). 
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Figure 7. CFIRB fluctuation amplitudes versus the 
CFIRB absolute level assumed. The hatched regions rep- 
resent the formal errors of the fits in Eq. ^| for 90 and 
170 /mi. A set of dashed lines marks <5B /B = 6, 7, 8, 10, 
12, 15 and 20%. 



At finer resolution (higher spatial frequencies) effects by 
individual galaxies have to be considered. An estimate 
of <5B /B (o r 5{vl v ) / (vl v )) can be performed following- 
Bond et al. ( 1986 ), assuming bia sed g alaxy clustering. 
Using eq. 6.35 from Bond et al. ( 1986|) , the fluctuation 
ratio is : 8(vl„) / (yl v ) = 0.05 x (l'/er) - 4 where a is the 
smoothing angle, and the beam profile is approximated by 
a Gaussian. Equations 6.24 and 7.2 in Bond et al. (1986) 
give a — 12" and 23" for 90 and 170 /im, respectively, for 
ISOPHOT. From this results 5{vl u )/{yl u ) = 10% at 90 /mi 
and 8(vl v )/(vl v ) = 7% at 170 /mi. 

At 170 /mi <5B /B = 7% 

B = O.SiO^MJysr- 1 equivalent 



provides 
to 



At 90 /mi the dependence of (5Bo 
even weaker than at 170 /im, therefore 
fined value of O.lliO.OSMJysr" 1 can be 



on Bo is 
a well de- 
determined, 
equivalent to 



The B value of l.l±0.3MJysr- 
v\ v = 37±10nWm~ 2 sr~ 1 ), is the same as the average sur- 
face brightness of the faintest 90 /im maps in our data set. 
The expected cirrus contribution in these fields - predicted 
from the ^0.8MJysr -1 value at 170 /mi by using a typical 
cirrus spectrum (modified black-body SED, v 2 emissivity 
law, 20K temperature) - is about 0.3MJysr _1 . This re- 



sult suggests that the l.lMJysr -1 is an upper limit and 
the most likely value of the CFIRB lies at the lower end 
of the uncertainty range, i.e. ~0.8MJysr _1 equivalent to 
i/I„ = 30nWm- 2 sr- 1 . 

Our results satisfy the three main criter ia of a CFIRB 
detection as proposed by Hauser & Dwek ( 2001 ): (1) We 
achieved a detection of a positive value of the fluctuation 
amplitudes (and therefore the absolute value as well) at 
90 and 170 /im at the 4<r level. (2) We removed the contri- 
bution by all known foreground (noise) emitters (zodiacal 
light, instrument noise, galactic cirrus). The properties of 
the spatial distribution related to the remaining noise term 
agree with that expected for the CFIRB fluctuations. (3) 
Due to the statistical nature of our method, the detection 
of the same well defined positive constant term in different 
areas of the far-infrared sky indicates the isotropy of this 
component. 

4.2.2. Comparison with other results 

The 170 /xm confusion noise Co = 15±4mJy is in good 
agreement with the 13mJy obtained by Juvela et al. 
(see the overview of Lemke et al. (2000)), and some- 
what lower than the values of 18mJy determined by La- 
gache & Puget (2000) from the analysis of the Maranol 



field and the 22mJy derived by Matsuhara et al. (2000) 



from the investigation of the Lockman Hole region. Our 
derived CFIRB value of 14±3nWm _2 sr _1 is the same 
as predicted from models by Pei et al. ( 1999| ). Also the 
COBE measurements calibrated on the FIRAS photomet- 



ric scale (Hauser et al. 1998, for an overview see Hauser & 



Dwek 2001) give quite similar levels (15±6 nWm sr at 
140/im and 13±2 nWm _2 sr _1 at 240/mi) despite the big 
difference in spatial resolution. 

The 90 /im confusion noise Co = 7±2 mJy is close (al- 
though somewhat lower) to the llmJy found by Mat- 
suhara et al. (2000). With regard to the 90 /im ab- 



solute value Schlegel et al. (1998) provided the same 
l.lMJysr -1 value as an upper limit for the CFIRB at 
100 /im, and also Lagache et al. ( 2000| ) and Finkbeiner et 
al. Q2000| ) reported a detection of a 100 /im CFIRB level 
of 23±6nWm~ 2 sr~ 1 and 25±8 nWm _2 sr _1 , respectively, 
which are close to our upper limit. 



4-3. Confusion limits for far-infrared space telescopes 

After the comparison of instrument and confusion noise 
values (see Fig. |^ and ||) it is obvious that even in the 
faintest sky regions the instrument noise is below the sky 
confusion noise by a factor of 2-3 for the C2_170 filter. 
This confirms the results of Herbstmeier et al. (1998), 



showing that the ISOPHOT C200 detector was limited 
by sky brightness fluctuations rather than by instrumen- 
tal effects. In the case of the Cl_90 filter the instrumental 
and confusion noise are similar in strength at the resolu- 
tion limit. Since our results confirm the applicability of the 
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H&B formula, the predictions made for other space mis- 
sions with telescopes of a similar size (SIRTF, Astro-F) on 
the basis of this formula can be trusted. Predictions made 
for the HERSCHEL satellite (3.6m mirror), however, still 
have to partially rely on assumptions, since we were not 
able to test the dependence of the confusion noise on the 
diameter of the telescope primary mirror. It remains to be 
demonstrated that the spatial structure of cirrus does not 
change below the resolution limit of ISO. 



5. Summary 

Using measurements obtained with ISOPHOT we investi- 
gated the properties of the sky confusion noise on a large 
sample of maps in the 90 < A < 200 /im wavelength range. 
We described the dependence of the sky confusion noise on 
the surface brightness for four selected ISOPHOT filters. 
We verified that the confusion noise scales as Ncx (B) 1-5 
for the resolution limit, independent of the wavelength. 
We also confirmed that, due to the dependence on the 
resolution parameter, N oc A 25 at A > 100 /im. Accord- 
ing to our results for cirrus fields with (B) > 5MJysr _1 
the Helou & Beichman ( 1990| ) formula predicts confusion 
noise values within a factor of 2. The scaling of the noise 
value at different separations between target and reference 
positions was investigated for the first time, providing a 
useful formula to estimate the confusion for different sepa- 
rations, too. At 90 and 170 /im a noise term with a Poisso- 
nian spatial distribution was detected in the faintest fields 
((B) < 3 — 5MJysr _1 ), which we interpret as fluctuations 
in the Cosmic Far-Infared Background (CFIRB). With a 
ratio of the fluctuation amplitude to the absolute level of 
10% and 7% at 90 and 170 /im, respecively, we determined 
the fluctuation amplitudes and the surface brightness of 
the CFIRB. The fluctuation amplitudes are 7±2mJy and 
15±4mJy at 90 and 170 /im, respectively. We obtained 
a CFIRB surface brightness of B = 0.8±0.2 MJysr" 1 
(y\ v = 14±3nWm~ 2 sr _1 ) at 170 /im and an upper limit 
of l.lMJysr" 1 (vl v = 37nWn:r 2 sr- 1 ) at 90/im. 
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